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Abstract. Symmetric logarithmic derivative (SLD) is a key quantity to obtain 
quantum Fisher information (QFI) and to construct the corresponding optimal 
measurements. Here we develop a method to calculate the SLD and QFI via anti¬ 
commutators. This method is originated from the Lyapunov representation and would 
be very useful for cases that the anti-commutators among the state and its partial 
derivative exhibits periodic properties. As an application, we discuss a class of states, 
whose squares linearly depend on the states themselves, and give the corresponding 
analytical expressions of SLD and QFI. A noisy scenario of this class of states is also 
considered and discussed. Finally, we readily apply the method to the block-diagonal 
states and the multi-parameter estimation scenarios. 


PACS numbers: 03.67.-a, 03.65.Ta, 06.20.-f 
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1. Introduction 

Quantum metrology has been going through a great development in recent years [1-16]. 
Quantum Fisher information (QFI) is a crucial concept in quantum metrology as it 
depicts the lower bound on the variance of an unbiased estimator for the parameter 
under estimation, according to the quantum Cramer-Rao theorem [17,18]. The definition 
of QFI is F := (L 2 ) = Tr(pL 2 ) [17,18], where L is so-called symmetric logarithmic 
derivative (SLD). Denoting the parameter under estimation as 6 , the SLD operator is 
determined by the equation 



( 1 ) 


Taking the trace on both sides of this equation, one can see that ( L) = 0. Therefore, the 
QFI is actually the variance of SLD operator, i.e., F = (A 2 L), with A 2 L := (L — (L)) 2 . 

The SLD operator has been studied for years [19-28]. It is important for two 
reasons. First, it is obvious that the QFI can be directly obtained when the SLD operator 
is known. Second, the achievement of quantum Cramer-Rao bound strongly depends 
on the measurement, namely, it can only be achieved for some optimal measurements. 
The eigenvectors of SLD operator are such theoretical optimal measurements [19-21] if 
they are locally independent of the parameter. Thus, the study of SLD operator could 
help us to construct or find optimal measurements for the achievement of the highest 
precision. 

The traditional method for the calculation of SLD operator is to expand it in 
the eigenspace of density matrix. We now denote the spectral decomposition of p as 
'^2fL 1 Pi\'f’i){f’i\, with pi, IVQ the ith eigenvalue and eigenstate of p, respectively. M is 
the dimension of p’s support. When p is positive definite (or full rank), M equals to 
the state’s dimension d. Utilizing the spectral decomposition form, the element of SLD 
operator can be expressed by [22-24,27] 



( 2 ) 


for i,j G [1, M\. Lij can be an arbitrary number when either i and j is larger than M. 
This method to calculate SLD operator is useful when the spectral decomposition of p 
is not difficult to obtain. However, for many cases, for example when the space of the 
density matrix is infinitely large, the diagonalization could be very tricky. 

Lyapunov representation is another method to obtain the SLD operator and applied 
in many scenarios [25-27]. The definition equation (1) is actually a special form of 
Lyapunov equation, indicating that SLD operator is a corresponding solution. In this 
representation, the SLD operator is expressed by [27] 



( 3 ) 


One advantage of this representation is that it is basis-independent. Generally, this 
representation is no easier to calculate than the traditional method. However, similar to 
the traditional one, Lyapunov representation would be very useful for some scenarios. 
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In this paper, we first review the Lyapunov representation and figure out that Eq. (3) 
is available for both full and non-full rank density matrices. Then we provide a further 
basis-independent expression of SLD operator based on the Lyapunov form. The new 
expression would be extremely useful when the anti-commutator between the density 
matrix and its partial derivative exhibits periodic properties. 

To show the advantage of this method, we apply it in a class of states showing a 
linear relation with their squares. This class includes all pure and all two-level states. 
We provide simple expressions of SLD operator and corresponding QFI via the given 
method for this class. Noise from the environment are widely exist in reality. The 
scenario for these states under depolarizing channel are considered. Moreover, we also 
discuss the block diagonal states and the multiparameter estimations. 


2. Lyapunov representation 


As the beginning, we first review the derivation of Lyapunov representation of SLD 
operator. Mathematically, Eq. (1) is known as a special form of Lyapunov equation. To 
solve this equation, one can construct a function 


f(s) = e~ ps Le~ ps , (4) 

which satisfies /(0) = L. The partial derivative of f(s) on s is d s f(s ) = 

— 2e~ ps (dgp) e~ ps . Integrating both sides of this derivative equation, one can obtain 

/»oo 

/(oo) - /(0) = -2 / e~ ps (d e p)e~ ps ds. (5) 

Jo 

When p is full rank, e~ ps trends to zero for s —> oo, indicating that /(oo) = 0. Thus, 
the SLD operator can be directly written in the form of Eq. (3). However, when p is 
non-full rank, /(oo) will not vanish. Reminding that M and d are the dimensions of p’s 
support and p respectively, then the limitation of e~ ps equals to diag {0^,1 cL-m} when 
s trends to positive infinite. Here 0 m is the M-dimensional zero matrix and t,i-M is the 
(d — M)-dimensional identity matrix. Correspondingly, we manually separate the SLD 
operator into four blocks as 


L 


Am 
B f 

JD d-M,M 


Bd-M,M 

Cd-M 


( 6 ) 


where the Hermiticity of L is applied. Utilizing this form, and based on Eq. (4), one 
can see that /(oo) = diag {0 m, Cd-M} ■ 

Meanwhile, on the right side of Eq. (5), in the eigenbasis of p, e~ ps is block diagonal, 
and the element [dgp]ij vanishes for i > M and j > M [23]. Then after the integral, 
the elements of the entirety for i , j E [d — M,d] keep vanishing. This means the block 
Cd-M is not involved in Eq. (5), and cannot be solved by this equation. Thus, Cd-M is 
actually undefined here. However, since Cd-M is also not involved in the calculation of 
QFI [22,23], therefore, it will not bring indeterminacy on the final expression of QFI. 
Based on this reason, we can simply take Cd-M = 0 for convenience. In this way, 


Quantum Fisher information and symmetric logarithmic derivative via anti- commutator s4 

the SLD operators for both full and non-full rank density operators can be uniformly 
expressed in the form of Eq. (3). 


2.1. Expanded representation 


Defining the anti-commutator p° as p°(-) = {p, •} and noticing the fact that the improper 
integral J 0 °° can also be written as Hindoo fj, Eq. (3) can be further expressed in the 
form 


L 


—2 lim 


\-S 


,n+1 


s-kxi' (n + 1)! 
n =0 v ’ 


(p°) n dgp. 


(7) 


This is a further basis-independent form of SLD. This formula would be very 
useful when the anti-commutator among p and its partial derivative exhibits periodic 
properties. In some cases, dgp is the eigen-operator of p°, i.e., p°dgp = adgp with a a 
real number. When a > 0, the SLD operator reduces to a very simple form 

L = -dgp. (8) 

a 

The simplest case here is the pure states. For a pure state, it is easy to see that a = 1 
as p 2 = p. Thus, the SLD operator for pure states is L = 2 dgp. 

Moreover, because of the equality 


0>T%>=£O’"W>/>) (9) 

m=0 


where C™ = n\/[m\(n — m)!], the SLD operator in Eq. (7) can also be written in an 
equivalent form 


L = 


-2 lim V V 


s 


>n+1 


n =0 m =0 


[n 


1 )! 


C™p m (dgp) p" 


( 10 ) 


This form of SLD operator could be very useful when p m (dgp) p n m are easy to calculate. 
Especially, when p commutes with dgp , above formula reduces to [27] 


L — p 1 dgp = (dgp) p \ 


( 11 ) 


where the equality Y2m =o ^n = has been applied. The detailed calculation can be 
found in the appendix. Based on this equation, the corresponding QFI reads 

F-Trfp- 1 ^) 2 ]. (12) 


One vivid example here is the parameter only shows up in the eigenvalues of p. In this 
case, dgp and p -1 are both diagonal in the eigenspace of p. Thus, the SLD operator is 
in the form of Eq. (11). More specifically, in the eigenspace of p. it is L i3 = Sijdg In pi. 

Recently, Correa et al. [31] studied a quantum thermometer prototype and 
demonstrated that after a full thermalization, the optimal probe state to access the 
maximum QFI is an effective two-level system with a highly degenerate excited state, 
of which the Hamiltonian can be written as H = Ai(/3) + ^2(/3)|£d)(£d|- 

Ai and A 2 here are temperature dependent energy gaps and |£j) is the eigenstate of 
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Hamiltonian for any i. In this case, since |£j) is temperature independent, the thermal 
state p = exp (—/3H)/Z is commutative with its partial derivative on temperature. 
Thus, its SLD operator can be expressed in Eq. (11). Specifically, it is L — — (Ai + 
fddpXi) Yl!i=\ |&)(£i| - (X 2 + /3d/3X 2 )l£d)(€dl- Substituting the expressions of Ai and A 2 in 
this equation, the QFI in Ref. [31] can be reproduced. 


2.2. Unitary parametrization 


A unitary parametrization process contains a large category of realistic parametrization 
processes. Recently, an alternative representation of QFI for unitary parametrization 
processes has been discussed [36-38]. For a unitary parametrization process, the 
parametrized state p = U(9)p in W (9), with U{9) a 0-dependent unitary matrix. The 
initial state p in is independent of 9. A key quantity in the alternative representation is 
a Hermitian operator Ft = i{dgU')U. All the information of parametrization is involved 
in this basis-independent operator. For a unitary process, the QFI can be expressed 
by [27,37] 


F = Tr (ftniej) , 


(13) 


where L e g = U^LU is the effective SLD operator. In this scenario, it is easy to find that 
p m (dep) p n ~ m = iUp 1 ^ [Ft, Pm] p n ~ m Uf Based on Eq. (10), the effective SLD operator 
can be written as 


J eff 


00 n , \n+l 

= -i2 lim V V 

s->oo (n + 1)! 

n =0 m =0 v J 




(14) 


For a pure initial state, as Pm[FL, pi n \p- m = 0, L e ff reduces to the known form L e ff = 
i2[Ft, Pm] [37]. When the equation {Ft, pf n } = 2p m FLp m is satished, p m and [7^,pi n ] are 
commutative. Then the effective SLD operator is L eff = i (Ft — p\ n FLpf n l ) . 


3. Application 

Now we apply Eq. (7) into a class of density operators, which share a common feature 
as below 

P 2 = «p-/3, (15) 

where a and /3 are real numbers. I 11 the eigenbasis of p, above equation is equivalent to 
pf = api~/3 for any i, which gives the solution p* = [a± y/a 2 — 4/3]/2. If only one of the 
solutions is positive, the density matrix is trivially proportional to the identity matrix. 
Thus, we only consider the situation that both solutions are positive, i.e., a > 2/3 > 0. 
Moreover, it is worth to notice that a, fd in Eq. (15) can either depends on 9 or not. 
Several well-known states, including all pure and all two-level states, satisfy this relation. 
From Eq. (15), it is easy to see that V = Trp 2 = a — (3 is the purity of p and satisfies 
dr 1 < V < 1 with d the dimension of p. From Eq. (7), the SLD operator for this class 
of states can be expressed by 

L = 1 [2 d e p + (dtP) p- 1 - dgalu] , 


(16) 
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where Im = diagjlM, Od-jvi} is the identity matrix on the support of p. The detailed 
derivation of this equation can be found in the appendix. If p is full-rank, Im is the 
identity matrix in the total d- dimensional space, i.e., Id- When p is non-full rank, i.e., 
det p = 0, p is the inverse matrix of p on the support. Based on previous discussions, 
we know that the elements of SLD outside the support are undetermined and they will 
not affect the value of QFI. Thus, for both full and non-full rank cases, Im can be 
replaced by Id, which means Ea. (16) can be simplified into 

L = ^ [ 2d g p + (dg/3) p~ l - d e a] . (17) 

When a, /3 are both constant numbers independent of 0, the SLD operator reduces 
to L — 2dop/a. Alternatively, a can be a constant and /3 is dependent on 6. Under 
this situation, the SLD operator reduces to L — [2 dgp+ (dg/3)p^ 1 ] /a. A well-known 
case here is the two-level states, which can be expressed in the Bloch representation 
p = (1 2 + r ■ cr) /2. Here a = (<j x , cr y , ct 2 ) t is the vector of Pauli matrices and r is the 
Bloch vector. Utilizing this representation, one can immediately find that p satisfies 
Eq. (15) with a = 1 and /3 = (1 — V)/2. Here V is the purity. Then the SLD operator 
is 

L = 2d e p- l -{d e V)p-\ (18) 

This is the general basis-independent expression of SLD operator for any two-level state. 
For a mixed state, the inverse of p can be written as 

p _1 = Y^j5 a y pT(T y ( 19 ) 

Thus, for a mixed state, Eq. (18) can be rewritten into 

L = 2d e p - Y^p a yP T(T y ( 20 ) 


Substituting the Bloch representation of p into above equation, the entire Bloch 
representation of SLD operator in Ref. [29] can be reproduced. 

Equation (17) is the general expression of SLD operator for all states satisfying 
Eq. (15). Utilizing this formula, the corresponding basis-independent expression of QFI 
reads 


F = 


cr 


2ctTr ( d e p ) + ( dg/3 ) Tr p - (2 M - 1) ( d e a ) d e /3 


( 21 ) 


The advantage of this basis-independent expression shows in two aspects. First, during 
the specific calculation, one can choose a convenient basis, in which d$p or p~ l are easy 
to express or calculate. Second, utilizing a basis-independent formula, the effects of the 
density matrix on the QFI are more distinct. 

The fact that Harmonic mean is less than the Arithmetic mean implies the 
inequality Trp -1 > M 2 . Meanwhile, it is easy to find Tr(dgp) 2 > 0. Thus, the QFI 
in Eq. (21) is bounded by the inequality 


F > L [(MdePf - ( 2 M - 1) (, 9,a) 3,13} . 


( 22 ) 
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This lower bound only depends on the coefficients and can be used to roughly evaluate 
the QFI. 

For the cases that a and fd are both ^-independent, the QFI reduces to F — 
2Tr (dgp) 2 /a. When a is a constant and fd is dependent on 6 , the QFI is in the form 

F = ^ [2aTr (, %p) 2 + (d 0 fd) 2 Trp" 1 ] . (23) 

From this equation, the QFI for any two-level state can be immediately obtained, which 
coincides with the equations given in Ref. [30]. For any mixed two-level state, the QFI 
is bound by the inequality F q > ( dgV ) 2 . 


3.1. States under depolarizing channel 


Most quantum states has to face the disturbance from the environment in reality. 
Depolarizing channel is very usual in quantum processes. Several quantum metrological 
problems of states under depolarizing channel have been discussed recently [32-34]. In 
this channel, the final state pf can be expressed by [35] 

1 — rj 


Pf = V Pin + 


d 


(24) 


where p m is the initial state and rj is the reliability of the channel. 

Now consider the situation that the initial state p m satisfies Eq. (15), i.e., p 2 n = 
CKinPin—An- Under this situation, the final state also satisfies Eq. (15) with the coefficients 
a = pa in + 2(l —p)/d, and fd = 7] 2 fd in + r](l ~v) a m/d+ (1 — p) 2 /d 2 . In this way, the SLD 
operator for the final state can be directly obtained by substituting the specific formula 
of a and fd into Eq. (17). 

If ai n and An are both d-independent, so will a and fd, the SLD operator then reads 

2 dp 


La — 


-deph 


(25) 


dr/a m + 2(1-77) 

One example for this case is a degenerate mixed state, i.e., p in = Y/iiLi 
where N (N < d) is the degeneracy and ('ipi(9)\'if>j(9)) = 5ij. It is easy to see that 
Pin = Pin/ Ft, satisfying Eq. (15). Taking 6 as the parameter under estimation, the SLD 
operator can be directly written as 

dp 


La — 


dp + 2N(1 — rf) Arj' 




(26) 


where L ini = 2<9 e (|'A.)(A|) is the SLD operator for |A)(A|- For a pure initial state, i.e., 
N — 1, the SLD operator reduces to the known form discussed in Ref. [34], 

Alternatively, the reliability p can also be the parameter under estimation. In this 
case, it can be checked that pf commutes with <9 r ?Pf- Thus, for any form of input state 
Pin, the SLD operator here is always in the form of Eq. (11), namely, L v = p^ 1 c> 7 ? pf. 
Denoting the spectrol decomposition of p in as Aj|Aj)(Aj|, the SLD operator can be 
expressed by 

L ” = A + 


( 27 ) 
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From above equation, one can see that the eigenbasis of initial state could be the optimal 
measurements to access the quantum Fisher information. The corresponding QFI is in 
then in the form 

(dXi - l ) 2 


f = 

v d[q(d\i - 1 ) + 1 ] 

If the initial state is a pure state, i.e., |A)(A|, the SLD reduces to 


L d 
11 (1 - rj)(r]d+ 1 - rj) 


|A)(A| 


1 — rj 


Id, 


and the QFI is 


F _ ("-D 2 

1 (drj + 1 — if) 2 


(28) 


(29) 


(30) 


From this equation, it can be found that with the increase of r), the value of F, n 
monotonously reduces. 


3.2. Block diagonal states 

The block diagonal states are widely used and discussed in quantum mechanics. One 
vivid example is the optical systems taking into account the superselection rules [39]. 
Generally, a block diagonal state can be written as p = ® [Li Pi- Here ® represents 
the direct sum. One can check that an available form of SLD operator here is block 
diagonal, i.e., L = ®" =1 L t , where L* is the corresponding SLD operator for p t . Consider 
the scenario that each block satisfies the equation p 2 = ctiPi — fa. It should be noticed 
that Trpj < 1 and the purity for p is V = JA cqTrp, — /%. In this scenario, each Lj 
satisfies Eq. (17). Thus, the entire SLD operator can be expressed by 

n i 

L = (J) — [2 d e pi + (deft) pf 1 - d e af\ . (31) 

i=l 

Moreover, if pi is a 2-dimensional block, it can be expanded via the Pauli matrices 
into pi = pit + Vi ■ a, where /q = Trpj/2, and r t is the Bloch vector for ith block. Then 
it can be found that pi satisfies Eq. (15) with a* = 2 pi and Pi = 2 p 2 — VJ 2. Here 
Vi = Trp 2 . With these coefficients, the SLD operator for ith block reads 

Li = — ( d g pi + fipf 1 ~ d e pi) , (32) 

where the coefficient Q = 2 PidePi — dgVi/A. If det pi = 0, Q vanishes. All the X states 
can fit in this scenario. 


3.3. Multiparameter estimation 

In multiparameter estimation, the quantum Fisher information matrix V is also defined 
via the SLD operators, i.e., 

•L ij ■ 2 Lqj}), 


(33) 
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where Lg i(j) is the SLD operator for parameter O^y For any state satisfying Eq. (15), 
Tij can be expressed by 


■Fij 


- oTr { dtp , djpj+di/3 ( dj/3 ) Trp 1 



(i diadj/3 + 


djadi(3 ) 


.(34) 


Here dyy represents the partial derivative on dy 3 y Obviously, the diagonal element of 
T reduces to the form in Eq. (21). For the cases that a is constant and fd is dependent 
on the parameters, T V] is in the form 


Fij = [aTr {dip, d 3 p} + d l d (djfd) Tr p . (35) 

Especially, for a two-level state, this equation reduces to 

B VB V 

J- q , y = Tr({^,^}) + ^-. (36) 


This is the general basis-independent expression of quantum Fisher information matrix 
for any two-level state. 


4. Conclusion 

In summary, we first reviewed the Lyapunov representation of the SLD operator and 
showed that this representation is available for both full and non-full rank density 
matrices. Furthermore, based on the Lyapunov representation, we gave a further method 
for the calculation of SLD operator. This method is particularly useful for those states 
whose anti-commutators with their partial derivatives exhibits periodic properties. 

As an application of the given method, we discussed a class of states, which have 
a linear relation with their squares. The corresponding analytical expressions of the 
SLD operator and QFI are provided via the method. Furthermore, we discussed the 
depolarizing channel scenario of these states and extend our discussion to the block 
diagonal states. For multiparameter estimation, the quantum Fisher information matrix 
is also analytically given. 

The calculation of SLD operator is an important topic in theoretical quantum 
metrology. We hope this work may draw attention in the community to studying more 
methods to obtain the SLD operators for various scenarios. 
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Appendix A. SLD for states commuting with their partial derivative 


In the following we give the detailed calculation of the SLD operator for the states 
commuting with their partial derivative. From the equation 

\ 71+1 


L = -2 lim V V 

s-s-no ' 


\7 

S) 


s-kxi ^^' (n + 1) 

n= 0 m= 0 v 1 


-Cffp m {dgp) p n 


one can see that when p commutes with dgp, this equation can be rewritten into 

00 ('_2s') n+1 

L — — lim Y A- > -—p n {dgp). 

(n + 1)! 

n =0 V ’ 

Remind that the spectral decomposition of density matrix is in the form 

M 

p=ypi\A){Ai 
2=1 


(A.l) 


(A.2) 


(A.3) 


where Pi and | iff) are ith eigenvalue and eigenstate of p, respectively. M is the dimension 
of p’s support. In this representation, the SLD operator is 

A (—2s) n+1 M 

lim > - 

s—^oo ■ 


L = 


—' (n + 1)! *— 

n =0 v ’ i=l 


y P m)(Hd s p 


M 


1 


y — lim (e 2spi - l) \A){A\deP 

L -* 71- O 


M 


2=1 

1 


Pi 


y~\ ^i)(A\dep- 

Pi 


(A.4) 


2=1 
yM -1 


It is known that Yli=iPi l'0«)('0i| is dehned as the inverse matrix of p on the support. 
Thus, the SLD operator can be finally expressed by 

L = p^dgp, (A.5) 

where p~ x is the inverse matrix of p on the support. 


Appendix B. Detailed calculation for the application 

In the application, for any state satisfying the equation 

p 2 = ap - £, 

one can obtain the following relation 

p° {d g p) = ad g p + {d g a) p - d e /3. 

Based on this equation and Eq. (B.l), the nth order term is in the form 

72—1 

(p°) n dgp = a n dgp + {pdga - doff) a^ 1 y 

m =0 



(B.l) 

(B.2) 


(B.3) 
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Submitting this equation into the equation 

(~s 


L = — 2 lim 


s-K» ( n -|_ 1)1 
n=0 V ’ 


(p°) n dep 


and since a > 0, the SLD operator can be expressed by 
2 

L = —d 9 p — 2 ( pd e a — dg/3) x 
a 


oo / \n+l n —1 /o 

. 7 

s->oo-^ (n + l)! ^ \ a 

n=0 y ’ m =0 v 


lim ) 

.<?—Vnn • ^ 


Utilizing the spectral decomposition of the density matrix, the term 

°o / \ n+l n -1 /o 


M oo f \ n+l n—1 / ^ \ m 


iim y^ y 

M oo 

ii m y—-—y 


Cl 


-E 


m=0 
oo / \ n+l 


a 


#!+)(+! 


(y r 


s->oa — 2 p; ' (n + 1)! 

i=i y1 71=0 v ’ 

M 


a 


n— 1 


l - 


2p, 

a 


lim y- 

s-loo^-j' q; — 2p.j 


-(e- sa -l)—— (e“ 

a v ’ 2pi v 


1 ) 


!+)(+ 


I+)(+!, 


(B.4) 


(B.5) 


where the equality YYYm=o x ™ = (1 — x n )/(l — x) has been applied. Since all pt here are 
larger than zero, above limitation reduces to the form 

M 




2=1 


(B.6) 


where p 1 is the inverse matrix of p on the support. Finally, the SLD operator for any 
state satisfying Eq. (B.l) can be expressed by 


L = — [2 d e p + (d e /d) p 1 - d e aI M ] , 


a 


(B.7) 


where Im is the identity matrix in the support of p. 


Appendix B.l. Calculation of QFI 

Based on Eq. (B.7), the quantum Fisher information F = ( L 2 ) can be directly calculated 
1 


as 


F = 


a 2 L 


4((<%p) ) - 4 {pea) (d g p) + (dg/3) Tr (p *) 


- 2 M (dga) {de!3) + {dea)‘ 


(B-8) 


where ({dop, p 1 }) = 2Tr (dgp) = 0 has been applied. Denoting the purity of p as V, 
i.e., V = Trp 2 , one can see that (dgp) = dgV/2, and 


(( Sepf) = It v\(p°dep)dep] 


(B.9) 
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Substituting Eq. (B.2) into above equation, there is 

((depf) = ^«Tr ( d e pf + ^ (d g a) ( d e V ). 

Therefore, the quantum Fisher information can be expressed by 


F = — 

cr 


2oTr ( d e pf - ( d e a) (d e V) + (d e aY 


+ ( d e /3f Tr (p x ) - 2 M ( d e a) ( d e f3) 

Moreover, since V = a — /3 here, above equation can finally be written as 


F = 


cr 


2aTr ( d g pf + ( d e (3f Tr p~ 


- (2 M - 1) (dga) (dgjd) 


(B.10) 


(B.ll) 


(B.12) 


Appendix B.2. Calculation of quantum Fisher information matrix 

For the SLD operators in Eq. (B.7), the element of quantum Fisher information matrix 
is 


iFij— 2Tr (p{dip, djp }) — djad{P — diadjV + diadja 
+ (diff) (djp) Tr p^ 1 - M ( djadifd + diddj/3). 
Since the first term in above equation can be rewritten into 
2Tr (p {dip , djp}) 

= i (diadjV + djadiV) + aTr ({dip, djp}), 
and d{P = d^a — dpi, can then be simplified as 


Fij = — 
a 2 l 


«Tr {dtp, djp} + d t fi (djp Tr p 

1 


-i 


-I M - - ) (diadjP + djadi/3) 


When i = j, above equation reduces to Eq. (B.12). 


(B.13) 
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